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The interaction with Ught weakens the superconducting ground state in classical superconductors. The 
situation in cuprate superconductors is more complicated: illumination increases the charge carrier density, 
a photo-induced effect that persists below room temperature. Furthermore, systematic investigations in 
underdoped YBa2Cu306+x (YBCO) have shown an enhanced critical temperature T^. Until now, studies of 
photo-persistent conductivity (PPC) have been limited to investigations of structural and transport 
properties, as well as the onset of superconductivity. Here we show how changes in the magnetic screening 
profile of YBCO in the Meissner state due to PPC can be determined on a nanometer scale utilizing 
low-energy muons. The data obtained reveal a strongly increased superfluid density within the first few tens 
of nanometers from the sample surface. Our findings suggest a non-trivial modification of the near-surface 
band structure and give direct evidence that the superfluid density of YBCO can be controlled by light 
illumination. 



Tuning the superfluid density of cuprate systems without introducing additional disorder into the structure is 
challenging. Illumination with photons is one tool to realize this. Photo -persistent conductivity caused by 
visible light is widely known in semiconductors such as doping-modulated superlattices of hydrogenated 
amorphous silicon (a-Si:H)^ or in low doped n-type Ge^. In these systems the recombination is inhibited by the 
spatial separation of electrons and holes, leading to long recombination lifetimes. In conventional superconduc- 
tors, short-time illumination with photons having energies larger than the Cooper pair binding energy causes 
transient pair breaking, leading to the creation of quasiparticles^. In this case, illumination reduces the energy gap 
and therefore weakens superconductivity. In contrast, in cuprates it has been demonstrated that mid-infrared 
femtosecond pulses induce transient superconductivity^. 

Remarkably, studies on YBa2Cu306+x with reduced oxygen content {x < 1) show a persistent increase in the 
electrical conductivity after continuous illumination with visible light, in thin films^'^ as well as in single crystals^. 
The critical temperature is shifted to higher values for a broad range of oxygen concentrations (0.4 < x < 0.9) 
up to AT^max — 15 K^. The magnitude of the shift strongly depends on the oxygen content x, the temperature at 
which the illumination is performed, and the photon dose^. YBCO samples with x < 0.4, which are close to the 
metal-insulator transition, even become superconducting under illumination^. The photo-induced effects are 
only present for photon energies higher than 1.6 eV. The changes in resistivity and persist until the temper- 
ature is raised beyond 270 K. Accordingly, in the literature the effect is referred to as photo-persistent conduc- 
tivity. Above 270 K the resistivity relaxes back on time scales of hours to days, depending on the temperature^. 
This time scale is similar to that required for YBCO samples to stabilize their properties at room temperature after 
deoxygenation of the sample at high temperatures The PPC is fully reversible and reproducible. If infrared (IR) 
light is used for illumination an increase in conductivity is observed as well, but with a much smaller magnitude^ \ 
But if IR illumination is performed after illumination with visible light, the PPC is reduced by about 2-10% even 
at low temperatures. 

Different origins have been proposed to explain the PPC, which can be generally divided into two main lines: 
(a) Photo-induced charge transfer and (b) photo-assisted oxygen ordering. The underlying mechanisms can only 
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be understood by taking into account the anisotropic structure of 
YBCO. In contrast to most other cuprate systems, YBCO consists of 
CUO2 planes separated by CuO^ layers which form Cu-O chains 
along the crystallographic b axis at higher doping levels x (Fig. 1). 
At low X the oxygen atoms in these CuO^ layers are isolated from 
each other^. Longer Cu-O chains form at higher x, leading to the 
formation of localized holes in the CuO^ layers. This results in a 
charge transfer of electrons from the CUO2 planes to the Cu-O chains 
increasing the number of mobile holes in the CUO2 planes. Therefore, 
the ordering of oxygen in the CuO;^ layers plays a major role for 
conductivity as well as for superconductivity in YBCO. A longer 
average chain length leads to a higher charge carrier density in the 
CUO2 planes and to a higher T^}^ on the underdoped side of the phase 
diagram. 

The scenario discussed above is crucial to an understanding of the 
proposed PPC mechanisms. In the case of photo-induced charge 
transfer, the absorbed photons create electron-hole pairs. The elec- 
trons are trapped at defects such as oxygen vacancies, preventing 
recombination. The holes are transfered to the CUO2 planes where 
they contribute to the conductivity and, at low temperatures to the 
superconductivity. Thus mechanism (a) explains the increase in con- 
ductivity and the shift in T^. However, it does not provide a natural 
explanation for the effects of IR light^\ The photo-assisted oxygen 
ordering mechanism (b) assumes that the photons initiate a reorder- 
ing of oxygen in the CuO;^ layers, so that longer Cu-O chain frag- 
ments are present. This leads to an enhanced charge carrier transfer 
to the CUO2 planes where superconductivity takes place. The pro- 
posed mechanism agrees with the observation that the time depend- 
ence of the resistivity relaxation is similar for photo -induced changes 
and for oxygen reordering after annealing. However, mechanism (b) 
alone fails to explain the IR experiments and why the oxygen atoms 
move under illumination. Therefore, a combination of both mechan- 
isms was proposed^ \ 

In the combined model the trapped electrons at the oxygen vacan- 
cies act as local perturbations that modify the local electric field 
distribution. The induced dipole moments cause the movement of 
the oxygen in the CuO^^ layers, resulting in a lengthening of the chain 
fragments. A simplified theoretical modeP^ agrees with this inter- 
pretation. It confirms that local perturbations can affect the average 
chain length. Experimental^^ as well as theoretical^^ investigations 
have also shown that electrical fields influence the oxygen ordering 
and change T^. The combined mechanism accounts for the infrared 
quenching, predicts the correct time dependence of the photo-excita- 
tion process, and gives the correct scale of the changes in resistivity 
produced by illumination. 

Even though extensive structural studies were able to correlate the 
chain ordering to the charge transfer, leading to changes of the res- 
istivity and Tc, little is known how illumination affects the super- 
conducting state on the microscopic level. In our study we focused on 
a detailed investigation of the magnetic screening profile B(z) in the 




CuO^ layer 



Meissner state, utilizing low-energy muon spin rotation (LE- 
Results & Discussion 

LE-^SR is an unique and powerful technique to measure non-trivial 
B(z) on a nanometer scale in a wide variety of superconducting single 
crystals^^, thin films^^'^°, and heterostructures^^'^^, allowing also the 
observation of non-local effects^^"^^ Positively charged muons are 
slowed down and implanted into the samples at different mean 
implantation depths z (Figs. 2 and 3). The Larmor frequency of the 
muons, Wi^(z), directly related to the local magnetic field B at the 
muon stopping site (col = with the muon gyromagnetic moment 
= 2n ' 135.5 MHz/T), is measured via the decay positron. For a 
semi-infmite sample, the London equation 



d^Bjz) 
dz^ 



(1) 



yields a magnetic penetration profile B (z) = B^^^ • exp (—z/li), for 
the boundary condition B{z = 0) = ^ext- Therefore, measuring B {z) 
allows one to determine the magnetic penetration depth Xj^ and 
thereby the superfluid density Usccl/Xl The process cycle and 
the conditions of the LE-^SR measurements are presented schem- 
atically in Fig. 3. The in-situ illumination setup is described in detail 
in the supplementary section SI. 

Two sets of YBCO thin films (x = 0.42, 0.67) with twinned crys- 
tallographic a and b axes were investigated along with a detwinned 
single -crystal mosaic (x = 0.67). The YBCO thin film sets show shifts 
of ATc = 6.4(6) K{x = 0.42) and AT^ = 0.9(1) K{x = 0.67) after 
illumination with about 7 * 10^^ photons/cm^, in agreement with 
previous results^ (see Fig. 4). The magnetic screening profiles B{z) 
determined by LE-^SR show for the initial state (before illumination) 
an exponential decay for the single crystals and a cosh behaviour in 
the thin films (resulting from the exponential penetration of the 
magnetic field from both interfaces) as depicted in Fig. 2. After 
illumination, all magnetic screening profiles are shifted to lower 
values and exhibit a pronounced change in shape. The corresponding 
superfluid density Us is enhanced in all sample sets investigated. 

Two different scenarios for the PPC change of as function of the 
depth z can be envisaged: (i) The additional charge carries are mobile 
along the lattice c direction, hence they are distributed homoge- 
neously over the whole sample. The superfluid density is raised, 
but is still constant as a function of depth z. Thus only the absolute 
value of is changed. Since the additional superfluid density is 
spatially averaged over the volume in scenario (i), the PPC changes 
of would be much less pronounced or more likely to be absent in 
single crystals compared to thin films. However, this does not 
account for our results, since a change in B{z) is observed in the 
YBCO ortho-VIII single crystals (see Fig. 2c). (ii) The PPC charge 
carriers are pinned to the layer where they are created. Here, the 
changes in the superfluid density would be proportional to the light 
penetration profile: Us increases mainly close to the vacuum interface, 
since the light penetrates exponentially into the sample. A depth- 
dependent superfluid density leads in turn to a modification of the 
Meissner screening field profile B(z). Our results agree with the 
second scenario. At small depths z the applied magnetic field is more 
strongly shielded compared to the initial state. At higher z values the 
B(z) curves are almost parallel. We therefore describe the penetration 
profiles of the illuminated samples by: 

' '-^^^Uz), (2) 



dz2 



/It Xt, 



Figure 1 | The crystal structure of YBa2Cu306+j£ 



where l/^^ is proportional to the initial superfluid density before 
illumination ns,o, l/ylp corresponds to the photo-induced superfluid 
density at the vacuum interface /is,iiium» and C is the characteristic 
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Figure 2 | Modifications of the magnetic screening profile due to illumination. The magnetic penetration profile B{z) along the crystallographic c 
axis normalized to the appUed magnetic field Bgxt of thin-film YBa2Cu306.42 (a), thin-film YBa2Cu306.67 (l>)> and detwinned YBa2Cu306.67 single crystals 
in the ortho-VIII phase for Bgxt applied parallel to the crystallographic a axis (c) and for 5ext parallel to the h axis (d). The applied magnetic field 5ext and 
the temperature are given in each figure. The presented field profiles correspond to the initial state (blue diamonds, before illumination), to the 
illuminated state (red circles, after illumination) and the recovered state (black triangles, after heating up to room temperature for about one day). The 
pronounced changes of the B{z) profiles (in magnitude and shape) due to visible light illumination are caused by a substantial increase of the superfluid 
density on the nanometer length scale to the vacuum interface. 



dph=1022 photons/cm^ 
I =80mW/cm2 
A = 400(10) nm 




270K 



Figure 3 | Schematic diagram of the experimental //SR process cycle. Left panel: After zero field cooling (ZFC) the sample an external magnetic 
field Bext/ A^o <^ci was applied within the a^-plane. The interface to the vacuum is labeled with z = 0. The screened magnetic field values B at different 
depths z along the crystallographic c axis were determined via the Larmor frequency of the implanted positively charged muons (red circles) col = y^B at 
10 K for the YBCO thin films and at 5 K for the YBCO ortho-VIII single crystals. The muon spin precesses at Wl around the local magnetic field B present 
in the sample. The Larmor frequency col is extracted from the muon-spin polarization function P(t) measured via the decay positron. Center panel: 
The LE-/zSR measurements were followed by the in-situ photo illumination of the samples for ~ 3 days (4-5 * 10^^ photons/cm^) at 270 K with a 
high-intensity LED source. Right panel: Afterwards the samples were ZFC to 10 K/5 K, repeating the measurement of the screened magnetic field profile 
as described above. To return the YBCO system to its initial state it was kept above room temperature for about one day. The resistance measurements 
have been performed under comparable conditions (see supplementary section SI). 
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Figure 4 | Shifts in the critical temperature AT^ from resistance measurements. Resistance R versus temperature T of YBa2Cu306.42 (a) and 
YBa2Cu306.67 (h) thin films for the initial state (blue diamonds, before illumination), the illuminated state (red circles, after illumination at 270 K), and 
the recovered state (black triangles, sample kept above room temperature for about one day). The critical temperature was determined as the 
temperature where R exceeds twice the RMS of the noise (dashed line) during heating of the samples. The shifts of the resistivity curves to lower values, 
after illumination at 270 K, reflect an increased charge carrier density at all temperatures T < 270 K. The arrows indicate the shifts in the transition 
temperatures ATc due to visible-light illumination. 



length over which the changes of take place. The solution of this 
phenomenological model, used for fitting, is described in detail in the 
supplementary section S3. The estimated values of C obtained from 
the fiSR data are of the order of the optical light penetration depth. 
Taking ellipsometry measurements of our thin film samples in the 
pristine state into account we find an average light penetration depth 
of C = 60(10) nm, consistent with previous studies on single crys- 
tals^^'^^, as presented in the supplementary section S2. Hence, we used 
this value for the detailed analysis of n^iz). 

The substantial changes of due to illumination are more pro- 
nounced for strongly underdoped samples as observed in the thin 
film sets: The measured change of the superfluid density at the 
vacuum interface /is,iiium/^s,o = 1 + {^J^pY is much larger for x = 
0.42 [n3,iiium/^s,o = 3.9(2)] than for x = 0.67 [^3^^5,0 = 1.45(3)]. 
In all sample sets investigated the PPC is fully reversible. After the 
thermal reset, the original Meissner screening profile is recovered as 
depicted in Fig. 2. 

Studies of single crystals were used to determine the role of disorder 
on the influence of the photo -induced changes. The investigated 
YBCO ortho-VIII single crystals consist of alternating Cu-O chains 
which are completely empty or full, forming a superstructure in the 
CuOx layers^^. Therefore, the amount of disorder is greatly reduced 
compared to our thin-film sample sets with equal doping level {x = 
0.67). Moreover, the single crystals are detwinned which allows us to 
investigate the anisotropy of the photo -induced effect. Assuming 
pinned charge carriers along the c axis [scenario (ii)], photo-persistent 
effects comparable to those in the thin-film sets are expected. In the 
ortho-VIII single crystal mosaic, the superfluid density at the vacuum 
interface is increased [ns,iiium/^s,o(ortho-VIII)= 1.7(1)] if the magnetic 
field is applied along the a axis (shielding currents flow along the Cu- 
O chains). The observation agrees with our thin-film results, implying 
that the Cu-O chains are also lengthened in less disordered single 
crystals due to illumination or showing that the order is not as perfect 
in the near-surface region compared to the bulk. The slighfly larger 
effect in the single crystal mosaic compared to the thin-film set of 
equal doping could result from the 15% higher photon dose due to the 
longer illumination time (see supplementary table SI). When the 
magnetic field is applied parallel to the b axis no effect on B{z) is 
observed. Within the experimental uncertainty the change of the 
superfiuid density at the vacuum interface /is,mum/^s,o has to be smaller 
than 1.06 for 5ext| |^-axis. 

The anisotropic PPC behaviour can be related to the in-plane 
anisotropy present in YBCO. The in-plane London penetration 



depths were determined as described in Refs. 18, 19 and found to 
be 4 = 177(5) nm {B^^tW and ^ = 158(2) nm (5ext||a). The 
anisotropy between la and li, is well known from optimally doped 
YBCO''-''. The ratio IJh = 1.12(3) of the ortho-VIII single crystal 
set is only slightly smaller compared to the values of optimally doped 
YBCO determined by various methods (^SR: XJlh ^ 1.16-1.19)^^ 
although it strongly depends on the doping level and the amount of 
oxygen disorder. The large in-plane anisotropy is discussed in the 
context of multiband effects^ \ This proximity model takes into 
account a contribution of the Cu-O chains to the superfluid density 
by assuming a hybridization of plane and chain wave functions 
within a three-band tight-binding Hamiltonian. Disorder within 
the Cu-O chains localizes a small fraction of the chain electron states, 
yielding a higher >^^, but has no effect on 1/^. Since the illumination 
of YBCO rearranges the Cu-O chains the amount of disorder within 
the Cu-O chains is reduced, thus diminishing the fraction of localized 
chain electrons. According to the proximity model, the additional 
superfluid density only affects >^^, consistent with our experimental 
data (see Fig. 2c-d). 

The presented LE-^SR data show a significant change in the mag- 
netic penetration profile in the Meissner state of underdoped YBCO 
caused by an accumulated superfluid density Us close to the surface. 
The persistent photo-induced changes of n^iz) appear only on a 
length scale comparable to the light penetration depth. The illumina- 
tion stimulates a self-organization of the Cu-O chains in YBCO and 
hence reduces disorder. This may suggest that the reduced disorder 
leads to a non-trivial modification of the local band structure, prob- 
ably due to hybridization of CUO2 planes and chain layers, which in 
turn strengthens the superconducting ground state. Our results show 
the tremendous impact of nano-scale disorder on the superconduct- 
ing ground state of cuprate superconductors. 

Methods 

Sample preparation and characterization. The YBa2Cu306+x (YBCO) thin films 
were deposited on single crystal 1X1 cm^ (001) SrTiOs substrates by on-axis pulsed 
laser deposition from a stoichiometric target using a KrF excimer laser (/I = 248 nm). 
To grow YBCO layers with a thickness of 400(20) nm and 490(60) nm, 3000 and 4800 
laser pulses with a laser repetition rate of 5 Hz were used under an oxygen partial 
pressure of 0.3 mbar. After deposition at 800°C, the samples were cooled to 780°C, 
where the chamber was flooded with 400 mbar oxygen to achieve fully oxidized films 
with a = 90.0(5) K. More details on the PLD growth are described in Refs. 32, 33. 
The layer thicknesses were determined by measuring the height of a chemical etched 
edge, ranging from the surface to the substrate, using atomic force microscopy. To 
check the homogeneity of the samples the thickness was measured at several points. 

The reduction of the oxygen content was realized by a post-annealing process in a 
tubular furnace, described in detail in Ref. 34. In particular, the samples were annealed 
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Figure 5 | Muon stopping profiles in YBCO. The normalized stopping 
distribution n{z) for positively charged muons at different implantation 
energies for YBCO was simulated with TRIM.SP^*^. The lines are guides to 
the eye. 



independent background of uncorrelated events. More details on the LE-/iSR tech- 
nique are given in Refs. 16, 17. The experiments were performed in the Meissner state. 
After zero field cooling to 5 K, a magnetic field B^^Jiuq < H^i was applied parallel to 
the surface of the sample. The muon-spin polarization function P{t) is then given by 



P{t) = e-'''''/^ jn{z)cos[y^B{z)t + (^]dz, 



(4) 



where ^ is the initial angle of the muon spin direction relative to the positron detector. 
The gyromagnetic ratio of the muon is y^ = 2n ' 135.5 MHz/T. The depolarization 
rate o" is a measure of any inhomogeneous local magnetic field distribution at the 
stopping site. The muon stopping distributions n{z) (see Fig. 5) were simulated for 
energies in the range 5-25 keV using the Monte Carlo code TRIM.SP^*^. The reli- 
ability of these simulations has been studied in various thin films^^'^''. The mean 

stopping depth z = j n{z) z dz of the muons in YBCO is in the range 25 - 106 nm. 

The corresponding mean magnetic field is B{z)= j n(z) B{z) dz. To analyse the 

data, the penetration profiles given in supplementary section S3 were used for B (z). 
The penetration depths /II were determined from global fits with the fiSR data analysis 
software package musrfif'", where fiSR spectra for different energies are analysed 
simultaneously. 

See supplementary information for full details on the illumination conditions, the 
ellipsometry measurements, and the penetration profile model. 



at a defined time and temperature in a flowing O2/N2 gas mixture. The process was 
concluded with a rapid cooling step to temperatures below 100°C in order to freeze in 
the oxygen concentration of the thin films. Annealing of the 400 nm thick samples at 
450°C for 40 min in 2860 ppm O2 atmosphere reduced the of the YBCO to 66.7(8) 
K (corresponding to x = 0.67). The 490 nm thick films with = 21(2) K (x = 0.42) 
were annealed for 70 min at 500°C in a 215 ppm O2 atmosphere. Purely c-axis 
oriented growth of YBCO was verified by 0 - 20 measurements using a Bruker D8 
Advance X-ray powder diffractometer with a Co anode. The resistivity measurements 
were performed with a four-point measurement on a cold finger cryostat in ultra high 
vacuum (about 10"*^ mbar) under comparable conditions used for the /iSR 
experiments. 

The YBCO ortho-VIII single crystals were grown by the self-flux technique 
described in detail in Ref 35. For the growth, BaZrOs ceramic crucibles were fabri- 
cated. After the self-flux growth the YBCO crystals were annealed in flowing oxygen 
for six days at 648°C, to set the oxygen content to x = 0.67. Afterwards the crystals 
were sealed in a quartz tube and annealed at 570° C for four days to homogenize the 
oxygen. During detwinning, the sample is heated to 180 - 250°C and squeezed with a 
pressure of about 100 MPa along one of the a/b axes under flowing nitrogen gas, 
causing the twin domains to reorient so that the a-axis (which has a shorter lattice 
constant) is along the direction of the pressure. For the formation of the oxygen 
vacancy ordered superstructure ortho-VIII, the YBCO single crystals were low tem- 
perature annealed at 40 - 42°C for five days. The final YBCO ortho-VIII single 
crystals had an estimated critical temperature of Tc ~ 66.5(5) K. One year later a = 
67.8(5) K was determined with a Quantum Design MPMS SQUID magnetometer. 
The increase of of about 1 K with time is natural in these crystals and is even 
expected, as the chain oxygen atoms slowly become better ordered even at room 
temperature and some of the interstitial oxygen impurities enter the Cu-O chains. 

Low-energy muon spin rotation. The LE-/iSR experiments were performed at the 
/^E4 beamline at the Paul Scherrer Institut (PSI, Switzerland). For the thin-film 
measurements a mosaic of four 1X1 cm^ large samples was glued with silver paint 
onto a nickel-coated aluminum plate. For the YBCO ortho-VIII measurements, a 
mosaic of 13 aligned single crystals, with an average size of about 2.1 X 1.7 X 
0.2 mm^, was used. The measurements have been performed in ultra high vacuum 
(about 10"^ mbar) on a cold finger cryostat. 

The /iSR technique uses positively charged muons (/i"^) as a local magnetic probe. 
To investigate the magnetic penetration depth profile B{z), the muons are slowed 
down using a solid Ar/N2 moderator, so that they stop on a nm scale in condensed 
matter. At the /iE4 beamline the low energy muons are produced at a rate of about 
10^ 3-1 16 y^iiii ~ 100% spin polarization^*^. The muons are implanted into the sample 
where they thermalize within a few picoseconds without noticeable loss of polariza- 
tion. There they decay after a mean fi~^ lifetime of t^^ = 2.197 ^s into a positron (e"^) 
and two neutrinos (v^, Vg): 

The positrons are preferentially emitted along the fi^ spin direction at the time of 
decay due to the parity violation in the weak decay. Therefore, the detection of the 
time difference t= t^ — between the implantation time and the decay time of the 
fi^ allows one to determine the temporal evolution of the muon-spin polarization P{t) 
using the expression^^: 

N{t)=No e-V^"[l+AP(0]+iVBkg, (3) 

where N{t) is the number of decay positrons at time t. The scale of the counted 
positrons is given by Nq. A is the observable decay asymmetry and A^skg is a time- 
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